THIN DOMAINS WITH DOUBLY OSCILLATORY BOUNDARY t 
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Abstract. We consider a 2-dimensional thin domain with order of thickness e which presents oscillations 
of amplitude also e on both boundaries . top and bottom, but the period of the oscillations are of different 
order at the top and at the bottom. We study the behavior of the Laplace operator with Neumann boundary 
condition and obtain its asymptotic homogenized limit as e — > 0. We are interested in understanding how 
this different oscillatory behavior at the boundary, influences the limit problem. 



1. Introduction 

In this paper, we analyze the behavior of the solutions of the Laplace equation with homogeneous Neumann 
boundary conditions 

-Aw e + w e = f in R e 

dw& n w (L1) 
— — = on dR 

dN e 

with f e £ L 2 (R e ) and N e is the unit outward normal to dR e . The domain R e is a two dimensional thin 
domain which presents a highly oscillatory behavior at the boundary and it is given as the region between 
two oscillatory functions, that is, 

R e = Uxx,x 2 ) € M 2 | xi e (0, 1), -eh( Xl /e a ) < x 2 < eg{x 1 /e)\, with a > 1. (1.2) 

where g, h : K — > K are C 1 periodic functions with period L\ and L 2 respectively (see Figure [TJ. Moreover, 
there exist constants ho > and h\,go,gi > such that < ho < h(-) < hi, and < go < <?(•) < g\- 

Observe that both the amplitude and period of the oscillations at the upper boundary, given by eg(x/e) 
are of the same order as the thickness of the domain. But, for the lower boundary, which is given by eh(x/e a ), 
the amplitude is of the same order e, while the period is of the order of e a , which means that we have much 
more oscillations at the bottom than at the top boundary, order of the lower oscillations is large than the 
order of the amplitude and height of the thin domain R e with respect to the small parameter e. 



The existence and uniqueness of solutions for problem ( 1.1 1 for each e > 0, is guaranteed by Lax-Milgram 
Theorem. We will analyze the asymptotic behavior of the solutions as e — > 0. 

Since the domain is thin, R e C (0, 1) x (— e/i(-), eg(-)), approaching the interval (0, 1), it is reasonable to 
expect that the family of solutions will converge to a function of just one variable and that this function will 
satisfy certain elliptic equation in one dimension with some boundary conditions. As a matter of fact, if the 
function h e (-) is independent of e, say h e {-) = 0, the limit equation is given by 

J -qoWxx + w = f(x), x € (0, 1) 

\ w'(0) = w'(l) = (L3) 

where qo = rpvj f Y * |l — §^~(?/i> J/2) | dy\dy 2 , and X is a convenient auxiliary harmonic function defined in 

the representative basic cell Y* = {(1/1,1/2) £ K 2 | < j/i < /, < y% < G(yi)}. 

The purely periodic case can be addressed by somehow standard techniques in homogenization theory, as 
accomplished in pflU]. See [5115] for general references in homogenization and [7] for reticulated structures. 
Observe that in this case the extension operators are very important for the convergence proof. 
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Figure 1. Thin domain R e 



Moreover, if we assume that g e (-) is independent of e, say g e (-) = g(-), and ho = min xe K{/i(a;)} then the 
variational formulation of the limit problem is: 

{(g(x) + ha) w x {x) ip x {x) + p(x)u)(x) <p(x)}dx = J f{x)<pdx, V^G-H' 1 (0,1) (1.4) 

where p(x) = g(x) + J^ 2 h(s) ds, for all x £ (0, 1), and the function f e (x) = f^\x) $ satisfies 

that f e — x /, w-L 2 (0, 1). We refer to [3] for details. In this work, we want to analyze the case where the 
thin domain is a region between two functions with different order of the oscillations. 

Our case is a combination of these two cases since both g € and h e are present. And we want to understand 
the effect of both terms at the same time in the limit equation. Notice that the techniques used to solve 
each case separately are different so we will need to combine both techniques to get the limit problem in our 
case. The main difference of the present work in relation to previous existing work in the literature, see for 
instance 6, 8] 0] and references therein, is that we allow two different order of oscillations in the boundary 
of the thin domain. 

In Section 2 we state the notation and the problem that we will study. Furthermore, we are going to 
construct an extension operator that will be very important in the proof of the convergence result. Finally, 
we state the main convergence result. 

In Section 3 we rigorously prove the convergence result. In order to do so, we combine two different 
techniques: we use an extension operator in the upper boundary and we define suitable rectangles in the 
lower boundary to apply the estimates that we obtained in Lemma |3.1| . 



2. Notation and statement of main result 

To study the convergence of the solutions of ( |1.1[ ) we first perform the change of variables (x, y) — > (x, ey), 
which transforms the domain R e into the domain £l c 

& = {(x 1 ,x 2 ) e M 2 | *i € (0, 1), -h(x 1 /e a ) <x 2 < g( Xl /e)}. (2.1) 
Under this transformation, we obtain the equivalent linear elliptic problem 



d Xl 2 e 2 dx 2 2 
du e 1 du e 
dx~^ + 7 2 dx- 



u t = f- in n c 



(2.2) 



/'i : — — — — v\ = on Oil' 

2 



where / e £ L 2 (Q e ) satisfies ||/ e ||L 2 (a«) < C, for some C > independent of e, and v e = {v\,v 2 ) is the 
outward unit normal to dft e . Observe that f2 e is not a thin domain anymore but there appears a factor 1 /e 2 
in front of the derivative in the x 2 - Moreover, the domain has very wild oscillatory behavior at the top and 
bottom boundary. 

For the analysis we will construct an extension operator for functions defined in the set fi e , but which will 
extend the function only over the upper part of the boundary. Hence, let us consider the following open set: 

& = {(aji.afc) G M 2 | xi £ (0,1), -h(xi/e a ) < x 2 < gA. (2.3) 

Lemma 2.1. With the notation above, there exists an extension operator 

p e £ £(L p {n e ), L p (h £ ))nc{w 1 - P {n £ ), w 1 *^)) 
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such that for any (p G W 1 ' p (Q e ), 

\\ P M\lp(SI') ^ KW^W LP (Q') 

dP e <p 



dP e ip 



and 



< K 



Lp(Q') 



&X 2 

dip 



< K 

LP(O e ) 



dip 



LP((l') 



dx 2 
dip 



dxq 



Lp(Q*) 



Lp(Q') d 



(2.4) 

(2.5) 
(2.6) 



where 1 < p < oo, K a constant independent of e and 77(e) = swp xeI {\g / e (x)\}. 



Proof. The extension operator is constructed with a reflection procedure over the upper boundary, as in 

Now, we state the convergence result: 

Theorem 2.2. Assume that / e G L 2 (fl e ) satisfies \\f e \\L 2 (n') < C with C independent of the parameter e 

f e (xi,x 2 )dx 2 - 



id that there exists f G L 2 (0, 1) such that f e — >■ /, w — L 2 (0, 1), where f e (xi) = f 91 



h(x 1 /e a ) 



Let u e be the unique solution of (2.2). Then, there exists uq G /f 1 (0, 1) such that if P e is the extension 
operator constructed in Lemma 2.1 we have \\P e u e — uqW^^^ and Uq is the unique weak solution of 
the Neumann problem 



where Y* is the basic cell 

Y* = {(2/1,2/2) G M 2 : < yi < L x and - h < y 2 < .9(2/1)}. 
The homogenized constant coefficients are defined by 



(2.7) 



f 91 1 f f dX 1 If 1 

4 - J_ ho *w ds = rJ Y A 1 - Wl ^) d ^ p -L- 2 J 



h(s)ds — ho, 



(2.8) 



where X is the unique solution (up to constants) which is L\-periodic in the first variable, of the problem: 

-AX = in Y* 
dX 



ON 



dX 



= onB 2 

g'(yi) 



(2.9) 



071 E>l 



{ 9N ^l+g'{ Vl f 

Bq is the lateral part of the boundary, B\ is the upper boundary and B 2 is the lower boundary of dY* 



Remark 2.3. If the non homogeneous term f e (xi,x 2 ) is a fixed function depending only on the first variable 

Li 



that is, f e (xi,x 2 ) — f(xi), it is easy to see that f{x\) = ( jf- +p)f(xi) o,nd therefore, (2.7) is the variational 



version of 



q 



Wxx + w = f(x), X G (0, 1) 

tu'(O) = w'(l) = 



(2.10) 



Notice that in case h(-) = 0, then p = and \y*T7Zi = \Y*\ Jy*(^ ~ J1/7) = ^° an< ^ we recover ( 1.3 > . 
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3. Proof of the main result 
The variational formulation of (2.2) is: find u e <E iJ 1 (fi e ) such that 



(3.1) 



Taking = u e in expression (3.1 ) and using that ||/ e ||L 2 (f2 e ) < C, we easily obtain the a priori bounds 

<9it' 



l|w £ |U 2 (o«), 



dx i 



and - 

l 2 (s>) e 



<9u e 



<9x 2 



L 2 (fi e ) 



< C. 



(3.2) 



If we denote by the standard extension by zero and by x e the characteristic function of f2 c , we may 



write (3.1 ) as 



du e 


dtp 


1 du e dip "i 


f f du f - 


dp 


1 du £ 


dip 


dx\ 


dxi 


e 2 8x2 dx2 J 


Jnt tdxi 


dx\ 


e 2 8x2 


dx 2 



( 3 - 3 ) 

where we divide the domain f2 c in two parts: one of them, fif_, carries all the oscillations and the other f2o 
is a fixed domain, that is, 



n € _ = {(x 1 ,x 2 ) G M 2 | xx G (0, 1), -/i(xi/e Q ) < x 2 < -&o} 
fi = {(xi, x 2 ) G M 2 | xi G (0, 1), -/i < x 2 < si}. 



(3.4) 



Before we start with the proof of the main result, let us state some relevant estimates on the solutions of 
certain elliptics problems, posed in rectangles of the type 



Q £ = {(x,y) G M 2 | -e a < x < e a , < y < 1}, with a > 1. 
As a matter of fact, for uq(-) G i/ 1 (— e Q ,e Q ), we define the function u e (x,y) as the unique solution of 



(3.5) 



= in Q e , 
dx e 2 dy 

u(x, 0) = Uo(x), on r £ , 

du 



(3.6) 



dv 



0, on dQ e \ T £ 



where v is the outward unit normal to dQ e and r e = {(x, 0) G M 2 \ — e a < x < e a }. 
We have the following, 

Lemma 3.1. With the notation from above, if we denote by uq the average of uq in T e , that is uq = 



l re 



f_ a uq(x) dx then there exists a constant C, independent of e and uq, such that 



1 /•£ 



J-e 



\u t (x,y)-u \ 2 dxdy < Ce a 1 ||r*o|||^(_e-, e «-) 



and 



du e 


2 


1 


du e 


dx 


i 2 (Q.) 


f "2 

e z 


dy 



i 2 (Q«) 



duo 
dx 



(3.7) 



(3.8) 



L 2 (- £ °, e °) 



Proof. See [3] for details. □ 



Proof of Theorem 2.2 The idea is to pass to the limit in (3.3) constructing appropriate test functions. First, 
we study the limit of the different functions that form the integrands of (3.3). 
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(a). Limit in the extended functions. Using the a priori estimate (3.2) and the results from Lemma 
2.1 we obtain that P e u e \n € iif 1 (ilo) and we can extract a subsequence of {P e u e \n g } C ff 1 (ilo), denoted 
again by P e u e , such that 

P e u f --uo w-H 1 ^) 

dP^ 2 (3-9) 

• s — L (U ) 



dx 2 



as e — > for some u G iJ 1 (f^o)- 



A consequence of the limits (3.9) is that ^0(^1,^2) does not depend on the variable x 2 . Moreover, we have 



that the restriction of P t u e to the coordinate axis x\ converges to uq. That is, P e u e \r — > uq s — H s (T) for all 
s E [0, 1/2) where T = {(xx, 0) <E M 2 \ X\ G (0, 1)}. Consequently, we obtain ||P e u £ - UolU^r) -> as e ->• 0. 
In view of the above limit, one has the L 2 — convergence of P e u e to Uq, that is 



\P e u e - uo|| L 2 ( fj e ) ~> as 6 -> 0. 



(3.10) 



In fact, on the one hand we have 



||P e u e (sBi,0)-«o(a:i)||^ (fil) 



1 rai 

/ |F e M £ (xi, 0) - w (xi)| 2 dx 2 dx\ 

J-h(x 1 /e a ) 

< C(g,h) \\P e u e -u Q \\ L 2 {r) -> as 0. 



On the other hand, 



|P e « 6 (a;i,a; 2 )-P e u e (a:i,0)|| 



L 2 (n«) 



1 rs 



< 



1 r a 



fl(ll/E») V JO 



fe 2 



-(xi,s) 



ds \x2\dx2dx1 < C(h,g) 



h(x 1 /e c ") 

8P f it 



P e u e (xi,x 2 ) — P € u e (xi, 0)| 2 dx\dx2 
2 



< eC{h,g) -> as 0. 



L 2 (O c ) 



Finally 



||P e w £ - u o|| L 2 ( o«) < ||P £ u e (^i,a;2) - P e w £ (xi,0)|| i2(Se) + ||P e w £ (a;i, 0) - Uo^OH^^ -> 0, 
as e — > 0. 

(b). Limit in the tilde funct ions . 

From the a priori estimates (3.2) we know that there exists a function £* G L 2 (Hq), such that, up to 
subsequences 



-Cw-L 2 (n ) and — 

0x1 0x2 



0s-L 2 (O o ); ase^O. 



(3.11) 



(c). Limit of y 6 - 

Let x be the characteristic function of the representative cell Y * . We extend x periodically on the variable 
yi G K and denote this extension again by x- Clearly, by construction, x £ (xi, X2) — x( x i/ e i ^2), for (xi,x 2 ) G 

Consequently, by the Average Theorem and the Lebesgue's Dominated Convergence Theorem we obtain 



x ee ^°e w*-L°°(Q ), where 0(x 2 ) := ±- [ ' X {s,x 2 )ds Vx 2 G (-h , 9l ). 

L i Jo 



(3.12) 



(d) Test functions. 

In order to construct appropriate test functions that will allow us to pass the limit in the variational 



formulation (3.3), we are going to need to define a partition of the unit interval [0, 1] which is related to the 
function h € and which will allow us to analyze in detail the effect of the oscillations at the bottom in the 
limit equation. Hence, denote by N e the largest integer such that N e L 2 e a < 1, where L2 is the period of the 
function h. Observe that N e ~ L^e - ™. Let 

h„ f — min 



(~), n = l,2...,N e 



(3.13) 
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and 7„, e € [{ n ~ ^)L 2 e a ,nL 2 e a ] a point where the minimum (3.131 is attained, that is, h{^^) = h n e where 
7„ i£ does not need to be uniquely defined. By extension, let us denote by 70, e = and 7^+1, e = 1. 

Note that the set {7o,e,7i,e, 7iV e +i,e} defines a partition for the unit interval [0,1]. Moreover, due to 
that h(-) is L 2 — periodic we have that h n e = ho for 71 = 1,2 ... , N e . 

We define now the test functions as follows. With <f> € ff 1 (0, 1), we consider ip e S H 1 (ft e ) defined as 



ip e (xi,x 2 ) 



n = l,2. 



(3.14) 



where is the rectangle Qn = {{ x i, x 2) | 7n.e < £1 < Tn+i.ej — hi < x 2 < —ho} and the function is the 
solution of the problem 

d 2 X% 1 d 2 X% . 

0, in Ql 



dx{ 
dX. 



e 2 dx\ 



* = 0, on flf£\T£ 



> X*(aj 1 ,x 2 ) = 0Ol)> onT^ 
where is the base of the rectangle, that is, = {(xi, —ho) : 7„ !e < x\ < 7„+i ie }. 
From Lemma 13. II we have 

2 



9X£ 


2 


1 




dxf 


H 







<Ce a-l ,1^,12 



i 2 (7n,e,7n + l,e)' 



Furthermore, since 



(/(xi,^) - (f>{xi) = ip € {xi,x 2 ) - tp e (x-L,0) 



d>p e 
9x 2 



(xi, s) ds, 



we have by (3.14) and (3.16) that 



(3.15) 



(3.16) 



V e - 011x3 (n«) as e ^ 0. 



(e) Passing to the limit. 



(3.17) 



We can now pass to the limit in (3.3) by making use of test functions tp e defined above. For this, we study 



the convergence of each term in (3.3). 
• First integrand: 



( du e dtp* ldu e dtp e ) , , f \ 1 1 / \ 7 7 

l<^(^ + ?aa^(9a^I J y Xl ' X2 > ( t ) y Xl > dx i dx 2 as e ->• 0. 



(3.18) 



Thanks to the choice of the test function (3.14) and the convergence (3.11), we easily get (3.18). 
Second integrand: 



f r du e dip' 1 1 du e dip e 1 

/ la a a ?«Zl«Z2 -> 

JsV l- oxi 0x1 e A ax 2 ox 2 J 



as e — >• 0. 



(3.19) 



From the definition of ip e , the Cauchy-Schwarz inequality and the inequality (3.16) we have (3.19). 
• Third integrand: 



X e P e u e ip e dx\dx 2 







puo(xi) 4>(xi) dx\ + / 9(x 2 ) Uo{xi) 4>(xi) dx\dx 2 as e — » 



(3.20) 



where the constant p is given by p = jt Jo* h(s)ds — ho. 



For this, note that we can rewrite the integral of the left side of (3.20) as 

/ x e -P c u e ip e dx\dx 2 = x c {PeU e — uq) ip e dx\dx 2 
Jci' Jn* 

+ / x c "o (ip e - 0) dxidx 2 + x e uo(j>dxidx 2 . 
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From (3.10) and (3.17), we have that the first two terms in the right hand side above go to 0. 
X e uo 4> dx\dx2 — uq (f>dx\dx 2 + / x tu o4>dxidx2 



Moreover, since 



= / Uq 



4> (h y-^j ~ ho^J dxi + I x^u-o&dxidx 



we get (3.20) from the Average Theorem and (3.12) 



• Fourth integrand: 



f e ip e dx\dx 2 



f(x\) (f>(xi) dxi as e — > 0. 



(3.21) 



From (3.17) and the hypotheses of the theorem we have (3.21). 



Therefore, using the convergences (3.19 ), (3.18 1, (3.20 ) and (3.21 ), we obtain the following limit variational 
formulation: 

{C(xi,x 2 )4>'{x 1 ) +9(x 2 )u {x 1 )(j)(x 1 )}dx2dx 1 + / pu (xi) (t>(x{) dx\ 



f{x 1 )<j){x 1 )dx l ,\l(l)&H x {Q,l) 



(3.22) 



with p = Jq 2 h(s)ds — ho. 

At this point the question is how to relate uq to £* . In the following subsection we will show a equation 
for£*. 

(f) Relation between £* and u . 

Let us consider the following families of isomorphisms T£ : A k \— > Y given by 



~/ v ,x\ — ekL\ 
T k {Xi,X2) = ( ,x 2 ) 



(3.23) 



where 



A l = {{xi,x 2 ) G I ekL 1 <xi< eLi(k + 1), -ho < x 2 < gi} and Y = (0,Li) x (-h ,gi)- 

with k G N. We can considerer extension operators P G C(H 1 (Y*),H 1 (Y)) D C(L 2 (Y*) , L 2 (Y)) , the proof is 
done in [3]. Using these operators, the isomorphism (3.23) and the unique solution of the auxiliary problem 
(2.9) we define uo k in {xx,X2) G A\ by 

,x\ — eL\k 



uj k (x 1 ,x 2 ) = x 1 -e(PXoT£(x 1 ,x 2 j)=xi-e(PX(- 



-,x 2 ) 



Observe that for any (xi,X2) G f2+ there is k such that (xi,x 2 ) G A\. Therefore, the function w e (xi,X2) 
ojf,(xi,x 2 ) is well defined and ui e G H 1 ^ 6 ^). We introduce now the vector rf = (r)l,r) 2 ) defined by 



du! c 

rjl(xi,x 2 ) = —(x!,x 2 ), (xi,x 2 )eU e + 



(3.24) 



where = {(xi,x 2 ) G M 2 : < x\ < 1 and — ho < X2 < g(x\/e)}. 

Taking into account the definition of X if we consider a test function ip G _ff 1 (il f ^) with ip — in 
neighborhood of the lateral boundaries, we get 

dip , e 1 dip \ 



Tll dx~ l +lh e 2 dx 2 J 



dx\dx2 = 0. 



(3.25) 



Then, with the variational formulation (3.1) and the identity (3.25) we can write: 



f r du e dip 1 du e dip , _, , i , , f 
Jq- I dxi dxi <l z 0x2 0x2 J Jn- 



dip 



1 5^ 



Vl dx^ + n2 e^ 0x2 



dx\dx 2 



J. M. ARRIETA AND M. VILLANUEVA-PESQUEIRA 

= [ x e f e *pdxidx 2 ,V(p £ H 1 ^). (3.26) 
We would like to pass to the limit in this expression. For this, we will construct appropriate test functions, 



which used in the identity (3.26) allow us to pass to the limit in all the terms. 

q, w*-L°°(n Q ), 



0, a-L^Qo); 



(g) Limit of uf and rj\. From the definition of cj £ , we have 

dx 2 

where 

1 - ^-(s,x 2 ) )x(s,x 2 )ds 



xi, s - L (fi ); 



q(x 2 ) 



(3.27) 



1 

Li Jo 



dyi 



See [2] for more details, 
(h) Function test. 

Let c/> = <P(xi) £ C^°(0, 1). We introduce the test function 



tp e (xi,x 2 ) 



X*(xi,x 2 ), (xi,x 2 ) £ Ctl n Q e „, ?i = l,2, 

<p(xi)u e (xi,x 2 ), (2:1,2:2) £ 0, € + = Q , 



(3.28) 



where cu e is defined above and, as in (3.14), Q e n is the rectangle Q e n — {{x\,x 2 ) \ j n<6 < x% < y n +i,ei —hi < 
x 2 < —ho} and the function is the solution of the problem 

d 2 Xt 1 d 2 Xf 



ax. 

dN' 



e 2 dx\ 



0, in 



2 = 0, on flf£\r< 



(3.29) 



, X*(xx,X2) = 0(xi)w e (xi,-/i o ), on T e r 



where is the base of the rectangle, that is, = {(xi, —ho) : 7„ !£ < X\ < 7„ + i j£ }. 
Moreover, we define the function X e (xi,x 2 ) — X^(xi,x 2 ) as (xi,x 2 ) £ Q c n D fil. 

and using the the properties of u> e we have that the function X e is H 1 (£l'L) and satisfies 



From Lemma 



3.1 



the following estimate 



dX e 


2 


1 


8X e 


dx\ 


L 2 (Oi) 




dx 2 



L 2 (f21) 



5(^(xi)w e (xi,-/i )) 
dxi 



L 2 (0,1) 



where C denotes a constant independent of e. Now, we can argue as in (3.17) and we obtain 

ll^-^Pwlia^j-^Oaac-^O. 



(3.30) 



(3.31) 



where Pui e is the function defined on the set {(xi, x 2 ) £ 



z\ £ (0, 1), —hi < x 2 < gi} using a extension 



operator obtained by reflection in the negative vertical direction along the line x 2 = —ho. .Indeed, since 
^ e {xi,x 2 ) - 4>{xi)Pw € {xi,x 2 ) = ip e (xi,x 2 ) - (j>{xi)w € {xi, -x 2 - 2h ) 
= ip e (xi,x 2 ) - ^ e (xi,-x 2 - 2h ) = f* 2 X2 _ 2ho ^{xi,s)d s ior(xi,x 2 ) £ 01 



we have by (3.28) 



\\r-^ e \\ L2{hc) <C(g : h) 
(i) Passing to the limit. 



dx 2 



C{g,h) 



L 2 (Q') 



duj £ 
dx 2 



i 2 (n;) 



C(g,h) 



dX £ 



dx 2 



as e -> 0. 



L 2 (n'_ 



Now we pass to the limit in the equality (3.26) considering the test functions (p — ip e and ip = </m e 
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First integrand: 



r 8u e 8ip e 1 8u e dip* ) 

< q — q h ^ ^ — q — >dxidx 2 -> as e -> 0. 

[j< I. axi oxi e- 2 OT2 ax2 J 



(3.32) 



Taking account the definition of -0 C , the Cauchy-Schwarz inequality and the estimate (3.30) we 
obtain the convergence (3.32 1. 
Second integrand: 



r 9u e 9?/> £ 1 8u e dip e e 8(6u e ) e 1 8(6u € ) ■> ^ 
I9xi9xi e 2 9x2 9x2 1 9xi ^ 2 e 2 9x2 J 1 2 

f ( ^ 86 86 i , , 

-> / K - — xx-q- — u ^axidx 2 as e ->• 0. 

./fin L ^1 ^1 J 



(3.33) 



From the definitions of r\\ and ip e the second integrand reduces to 

Therefore, using convergences (3.10), (3.11) and (3.27), we have (|3.33|. 
Third integrand: 

X e P e u € i\f dx\dx2 — > I p Uq(xx) 6{x\) xi dx\ + / 0(x2) tto(xi) 6(x\) xi dx\dx 2 , as e — > 0. (3.34) 



52' 



Following along the lines of the proof of the convergence (3.20) we have this convergence. 
• Fourth integrand: 

/ e -0 C dx\dxi — > / f{x\) <p(xi) X\dx\ as e — > 0. (3.35) 
fi e ./o 



Using the same computations as those made to derive (3.21) we obtain (3.35) 



considering the test functions ip = ip e and ip — 6u e . More precisely, we have 
86 



Now, by the convergences shown in (3.32), (3.33), (3.34) and (3.35), we can pass to the limit in (3.26) 
ddering the test functii 

r 86 86 \ 

L 8x\ 8x\ > 

where p and q are given by 



.rj — (/— — //,) f dx\dx<i + / puQ,6x\dxi+ j 6uq 6x\ dx\dx 2 = j f6x\dx\ V<?!> € Co°(0, 1) 

(3.36) 



P 



J i Jo 



h(s)ds — ho, q(x 2 ) 



L 



l Jo 



1 - d ^ L {s,x 2 ) )x{s,x 2 )ds. 



dyi 



Taking the test function 6xi in (3.22) we obtain 
8 



9xi 



£*(x 1 ,x 2 ) — — (6xi)dx2dxi + / puo(xi) 6(xi) xi dxi + / #(x 2 ) u$(xi) 4>{xi) x\ dx\dx 



n 



f(x 1 )6(xi)x 1 dxi (3.37) 



Due to C^r(<f>Xi) = Cxi-§£- +C4>, w e can rewrite (3.36) as 



(f*7q — {(j>xi) — 6£* - q- — u \dxidx2 + / pu Q 6x\ dx\ + j 9u 6x 1 dx\dx2 
o L 9xx 9xi J Jo Jn 



f6x x dx x MS e C£°(0,1). 



Therefore, it follows from p737| ) and ( |3T38| ) that, for all G C£°(0, 1) 

= ^" ^* + q^^uojdxxdx 2 = J |^* - g^</>|dxidx 2 



(3.38) 



(3.39) 
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With the definition of q given by (2.8) and performing an iterated integration in (3.39) we obtain 

)dx 1 = 



(j>(x x ) 

So, the equation satisfied by £* is: 



t, {x 1 ,x 2 )dx 2 - q 

-h 



dx\ 



V0 € Cg°(0,l) 



4 (xi,X2)dx 2 = q-. 

-h„ 



dx\ 



The last step is placing this last equality in (|3.22|). We get 
rl (du dip \Y 



f ( du dip 
J Q I dxi dxi 



+ '—^ - u ip + u ipp >dx 



}dxi = J ftpdxi, V^Glf^O,!). 



(3.40) 



Hence uq is the unique solution of (3.40), and we obtain that any convergent subsequence of {u c } tends 
to this unique solution. This complete the proof of Theorem 2.2 □ 
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